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Abstract 
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1 Introduction 



Recently, considerable interest |], ||, ||, f|, ^ in the universal seesaw mass matrix 
model H has been revived as a unified mass matrix model of the quarks and leptons. 
Suggested by the seesaw mechanism for neutrinos 0, the model was proposed in 
order to understand the question why the masses of quarks (except for top quark) 
and charged leptons are so small compared with the electroweak scale (~ 10 2 
GeV). The model has hypothetical fermions Fi in addition to the conventional 
quarks and leptons fa (flavors / = u,d,v,e; family indices i = 1,2,3), and those 
are assigned to f L = (2,1), f R = (1,2), F L = (1,1) and F R = (1,1) of SU(2) L x 
SU(2)i?. The 6x6 mass matrix which is sandwiched between the fields (fp,Fp) 
and (f R , F R ) is given by 

M 6 * 6 = ( ° mL I , (1.1) 

m R Mp 

where mp and m R are universal for all fermion sectors (/ = u,d,i>,e) and only 
Mp have structures dependent on the flavors /. For < A^ <C As, where 
Ap = 0(mp), A R = 0(rrip) and As = O(Mp), the 3x3 mass matrix My for the 
fermions / is given by the well-known seesaw expression 

Mf ~ -m L Mp X m R . (1.2) 

However, after the observation || of the heavy top quark mass m t ~ Ap, the model, 
at one time, became embarrassed, because the observed fact m t ~ 0(mp) means 
0(Mp 1 m R ) ~ 1. This problem was recently solved by Fusaoka and the author[|T|, 
and later by Morozumi et al. 0. If we can built a model with detiW^ = for the 
up-quark sector (F = U), one of the fermion masses m{Ui) is zero [say, m(C/ 3 ) = 0], 
so that the seesaw mechanism does not work for the third family, i.e., the fermions 
(u3L,U3 R ) and {uzr^Uzl) acquire masses of 0{mp) and 0(mp), respectively. We 
identify (u 3 p, U 3r ) as the top quark (t L ,t R ). Thus, we can understand the question 
why only the top quark has a mass of the order of A L . Of course, we can successfully 
describe JU] the quark masses and mixings in terms of the charged lepton masses by 
assuming simple structures of mp, m R and Mp. The model also gives an interesting 
phenomenology for neutrinos 0. 

In spite of such phenomenological successes, there is a reluctance to recog- 
nize the model, because the model needs extra fermions F. In most unification 
models, there are no rooms for the fermions F. For example, it has been found 0] 
that when the gauge symmetries SU(2)ixSU(2) j RxU(l)y xSU(3) c are embedded 



into the Pati-Salam@ type unification SO(10) -> SU(2) L xSU(2) /? xSU(4) P5 , those 
gauge coupling constants are unified at \i — A x — 6x 10 17 GeV [SU(4) is broken into 
U(l) y xSU(3) c at n = A R ~ 5 x 10 12 GeV]. However, in the SO(10) model, there 
is no representation which offers suitable seats to the fermions F L m = (1, 1,4)l/r 
of SU(2) £ xSU(2)^xSU(4)p5. Whether we can built a unification model in which 
the fermions F are reasonably embedded will be a touchstone for the great future 
of the universal seesaw mass matrix model. 

For this problem, there is an idea fTo|- We can consider that the fermions 



Fr (= CF R ) together with the fermions fi belong to 16 of SO(10), and also F[ 
together with f R belong to 16 of another SO(10), i.e., 

+ *S) ~ (16, 1) , (/r + F£) ~ (1, 16) , (1.3) 

of SO(10) L xSO(10) fi . The symmetries are broken into SU(2) L xSU(2) R xU(l) y xSU(3) c 
at fi = As and the fermions F have the mass term F^MfFr. 

In order to examine the idea (1.3), in the present paper, we investigate the 
evolution of the gauge coupling constants on the basis of the SO(10)lxSO(10)r 
model and estimate the intermediate energy scales Ar and As together with the 
unification energy scale Ax- Of the numerical results, especially, we interest in 
the value of k = Ar/Al, because the value is closely related to the neutrino mass 
generation scenarios as we discuss in the next section. The evolutions of the gauge 
coupling constants under SO(10)iX SO(10)^ symmetries have already been done by 
Davidson, Wali and Cho ||10[| , but their symmetry breaking patterns are somewhat 



different from that in the present model. We will investigate the possible intermedi- 
ate energy scales under a constraint Ar/As — 0.02 Q which were derived from the 
observed ratio m t /m c in the new scenario of the universal seesaw model jl], where 
the masses m t and m c are given by m t ~ and m c ~ (Ar/ As)Al, respectively. 

In Sec. 3, we investigate the case of the symmetry breaking SO(10) L x 
SO(10) fl -> [SU(5) x U(1)'] L x [SU(5) x XJ(1)'] R . We will ruled out this case, 
because the results are inconsistent with the observed values of the gauge coupling 
constants at \i — m z . In Sec. 4, we investigate the case SO(10) L x SO(10) fl — > 
[SU(2) x SU(2)' x SU(4)]i x [SU(2) x SU(2)' x SU(4)] fl . We will conclude that the 
case is allowed for the intermediate energy scale Ar ~ (10 1 — 10 6 ) GeV if we accept a 
model with Axl ^ Axr, where Axl and Axr are the unification scales of SO(10) L 
and SO(10) fl , respectively. Finally, Sec. 5 will be devoted to the conclusions and 
remarks. 



2 Neutrino mass matrix 



In the universal seesaw mass matrix model, the most general form of the neutrino 
mass matrix which is sandwiched between (z/ L , T> C R , Nl, N r ) and (z/£, ur, iV£, Nr) t 
is given by 



M 



12x12 








IT 







m R 
j 



m 



R 



m 



m R 

Mr 

M T D 



m L \ 

IT 

m R 
M D 
M L j 



(2.1) 



under the broken SU(2) L x SU(2) R x U(1)y symmetries. Here, we have denoted 
the Majorana mass terms of the fermions F£ and F R as Mr and Ml, respec- 
tively, because the fermions Fl and Fr are members of (1, 16*) and (16*, 1) of 
SO(10) L x SO(10) fl , respectively. The mass terms / ' lTUlFr and FLtriRfR are gener- 
ated, for example, by the Higgs scalars (126, 1) and (1, 126*) of SO(10) L x 80(10)^, 
respectively, while the mass terms j ' L m! L Fl and F R m' R f R must be generated by 
Higgs scalars of the type (16, 16*) of SO(10) L x 80(10)^. Therefore, in the present 
model, we do not consider the terms m' L and m' R , i.e., we take m' L = m' R = 0. (For 
the special case with m' L ~ m L and m' R ~ m R , see Ref. |nj.) Hereafter, we assume 
m L <^ m R M F . 

Our interest is in a mass matrix for the left-handed neutrino states vl- By 
using the seesaw approximation for the matrix (2.1), we obtain the 6x6 mass 
matrix for approximate (z/£, vr) states 



M 



6x6 







III 




Mr M D 

Ml M L 







ml 




m L M 22 1 ml m L M 21 1 mR 
m T R M^ 2 l m T L m T R M^m R 



',2.2) 



where 



/ Mr 


M D > 


"'-( 


[Ml 


M L j 





M 21 l 



M X2 

M 22 l 



;2.3) 



Afu = Mr — M D Ml l Ml , 
M 22 — Ml — M T D M R X M D , 



(2.4) 



A 



Mia = =M%- M l Md X M r ■ 

According as the cases (a) M L , M R > M D , (b) M L ,M R ~ M D , and (c) M L ,M R <t: 
Mp, we obtain the following mass matrix for the approximate vl states, 
(a) The case M L , M R > M D 

From M u ~ M R , M 22 ~ M L and M 12 ~ -M L M^M R , we obtain 



u <;<> . -^LM^ml m L Ml x M T D M R x m R 

' m T R M R l M D Ml x m T L -m T R M R l m R ' ' ' ' ° 



so that we get the mass matrix for approximate i/£ states 

M(z/ L ) ~ -m L Ml x m T L , (2.6) 

because of (M 6x6 ) n , (M 6x6 ) 22 > (M 6x6 ) 12 . 
(b) The case M L} M R ~ M D 
We consider the case 



<M| M ; Md ^0. (2.7) 
1 Ml M L 1 



(The special case that the determinant is zero has been discussed in Ref. ||.) 
Since we consider the case m L <C m R , we can use the seesaw approximation for the 
expression (2.2), so that we obtain 

M(v L ) ~ -m L M 22 l m T L + m L M 2 ^m R (m\M^m R y l m T R M{ 2 l m T L 

= -m L (M 22 - M 21 x M u Mi 2 )ml = -Tn L M£ l ml , (2.8) 

where we have used the relation Ml = (M 22 — M 2 ~ 1 1 MnM 1 ~ 2 1 ) _1 in the inverse 
expression of (2.3). Thus, we obtain the expression (2.6) for the case (b), too. 
Note that the 3x3 mass matrix for approximate vl states is almost independent 
of the structures of M D and M R in spite of 0(M L ) ~ 0(M D ) ~ 0(M R ). 
(c) The case M L , Mr < M D 

From M n ~ -MdM^M^, M 22 ~ -MlM R l M D and M 12 ~ M^, we obtain 
the mass matrix 



ylixli ^ , m L MB 1 M R M^- 1 ml -m L M^m R 

1 -m R M T D - x ml m T R Ml- l M L M D x m R 1 



The mass matrix gives three light pseudo-Dirac neutrino states p2[ if± ~ (vn ± 
v? R )/V2 (i = e,fi,r), because (M 6x6 ) n , (M 6x6 ) 22 < (M 6x6 )i 2 . This case has been 
discussed by Bowes and Volkas[TB|]. The case is very attractive phenomenologically, 
because the maximal mixing state between and v^ R can give a natural expla- 
nation for the recent atmospheric neutrino data|H|. The mass matrix M(v± D ) in 
the limit of m(t / f_^' D ) = m(vfl D ) is approximately given by 

M(v p ± D ) ~ -m L M^m R . (2.10) 

First, we suppose the following symmetry breaking pattern (hereafter, we 
will refer to it as the case A): 

SO(10) L x SO(10) fl 

I At = A X io 
[SU(5) x U(1)'] L x [SU(5) x U(l)1fl 

I ^ = A N 

SU(5) L x SU(5) R 

[SU(3) x SU(2) x U(1)] L x [SU(3) x SU(2) x U(1)] B 

1 fi = A S 
SU(3) C x SU(2) L x SU(2) R x U(l) y 

1 A* = Afl 
SU(3) C x SU(2) L x U(l) y , 
1 A* = A L 

SU(3) C x U(l) em . (2.11) 

At the energy scale /i = A N , the gauge symmetries XJ(1)' L x XJ(1)' R are completely 
broken, so that the neutral leptons Nl and N R acquire Dirac and Majorana masses 
of the order of A^. At \i — As, the remaining fermions Fl and F R (except for 
Usl and Us R ) acquire masses of the order of As by Higgs bosons $ (as we discuss 
in the next section), and SU(3) L x SU(3) fl and U(1) L x U(l) fl are broken into 
SU(3) L+ij , = SU(3) C and U(l) L+i? = U(l) y , respectively. If this scenario A is true, 
the mass matrices M L , M R and M D are of the order of A N , so that we suppose that 
the order of the neutrino masses m(ui) are given by 

m[v t ) ~ A\jA L ~ (A L A s /A R A N )m(ei) , (2.12) 



from the result (2.8) in the case (b), the neutrino masses are suppressed by a factor 
(A L /A R ) (A S /A N ) compared with the charged lepton masses m(ej). 
Next, we can suppose another symmetry breaking (case B): 

SO(10) L x 80(10)* 

I (J> = A x 

[SU(2) x SU(2)' x SU(4)] L x [SU(2) x SU(2)' x SU(4)] fl 

I M = A S 

SU(2) L x SXJ(2) R x U(l) y x SU(3) C 
I H = A R 
SU(3) C x SU(2) L x U(l) y , 
I H = A L 

SU(3) C x U(l) em . (2.13) 
If this scenario B is true, since Ml ~ M R ~ Md ~ Mg, we suppose 

m(^) ~ A 2 L /A S ~ (A i /A«)m(e i ), (2.14) 

so that the neutrino masses m{y^) are suppressed by a factor A L / A R compared with 
the charged lepton masses m(ej). 

What is of the great interest is to estimate the possible values of such inter- 
mediate energy scales A*, A s , and so on. 

Although the Bowes- Volkas model fl3j is very interesting, the model cannot 
apply to the universal seesaw model based on the SO(10) L x SO(10) fl unification, 
because the case Ml, Mr -C Md is not likely in the SO(10) L x SO(10) fl model, and, 
if it is possible, the relation (2.10) leads to the wrong prediction m{Pi) ~ ra(ej) for 
M n - M s - Mp (F ^ N). 



3 Case of SO(10) SU(5) x U(l) 

In the present section, we investigate the case A with the symmetry breaking 
pattern (2.11). At the energy scale \i = As, the symmetries [SU(3) x SU(2) x 
U(1)] L x [SU(3) x SU(2) x U(1)] R are broken into SU(3) C x SU(2) L x U(l) y by the 
following Higgs scalars $y: 

$2/3 ~ (3*, 1; 3, l)y= 2 /3, 
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$4/3~(3,l;3*,l) y= 4 /3 , (3.1) 

$ 2 ~ (l,l;l,l)y =2 , 

of [SU(3) x SU(2)] L x [SU(3) x SU(2)] R , where SU(3) C = SU(3) L+K , U(l) y = 
U(l) i+fl , and Y = Y L = Y R . Our interest is in the region A L < /j, < A X5 . 
Hereafter, we call region A L < /j, < A R , A R < /j, < A s , and A s < pL < A X5 regions 
I, II, and III, respectively. 

The electric charge operator Q is given by 

Q = J 3 L + l -Y' (Region I) , (3.2) 

l -Y' = I? + \Y (Region II) , (3.3) 

\y = l -Y L + l -Y R (Region III) . (3.4) 

We denote the gauge coupling constants corresponding to the operators Q, Y', Y, 
Y L , Y R , I L , and I R as g ern = e, g[, g u g 1L , g 1R , g 2L , and g 2R , respectively. The 
boundary conditions for these gauge coupling constants at /i — A L , /i — A R , and 
/i — As are as follows: 

a^(A L ) = a£(A L ) + j^'f 1 ^) , (3.5) 
jja'f \A R ) = a^(A R ) + ^\A R ) , (3.6) 

and 

l*i\As) = l^l(As) + l^(A s ) , (3.7) 

respectively, correspondingly to Eqs. (3.2), (3.3) and (3.4), where «j = gl/tot 
and the normalizations of the U(l)y, U(l)y, U(l)y L and U(l)y R gauge coupling 
constants have been taken as they satisfy a[ = a 2 L = &3, Q-il = Q-2L = ®?>l, 
and aii R = a 2R = a^ R in the SU(5) grand-unification limit and ot\ = a$ = 
in the SU(4) unification limit [a^ = a 2 L = ol 2 r in the SO(10) unification limit], 
respectively. We also have the following boundary conditions at /i — As and 
li = A X5 : 

a^(A s ) = a-l(A s ) + a^(A s ) , (3.8) 
«r2( A x5L) = ^l(A X5L ) = a^l(A X5L ) , (3.9) 
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«ii(AxBfl) = ®2r( a x5r) = a 3 R(A X 5R) , (3.10) 

where, for convenience, we distinguish the unification scale of SU(5) L , A 5 xl, from 
that of SU(5) H , A X 5R- 

The evolutions of the gauge coupling constants gi at one-loop are given by 
the equations 

^oii(n) = -^-b^ifi) , (3.11) 

where t = In fi. Since the quantum numbers of the fermions / and F are assigned 
as those in Table [l], the coefficients hi are given in Table |[ In the model with 
detMj/ = 0, the heavy fermions Fl and Fr except for U^l and U^r are decoupled 
for fi < A s and the fermions u^r and U 3 l are decoupled for fi < A R . In Table 
D, we have also shown the values of 6j for the conventional case without such the 
constraint detMjj = in parentheses. 

By substituting o:^2(Ax5l) = oc^I(Ax5l) with the relations at one-loop 

a^(A X5L ) = a 2 l(A s ) + In , (3.12) 



o£(Axbl) = «3l(As) + b 1 "^- In ^ , (3.13) 

27T Ac 



we obtain 



^l(As) - a£(As) + (bg 1 - b&')±- In ^ = . (3.14) 

27T A 5 



Similarly, from the condition a 1L {AxhL} = «2l(Ax5l), we obtain 

a£(As) - ocTt(As) + (6ff - b%)±- In ^ = . (3.15) 

By eliminating ln(Ax5z,/As) from Eqs. (3.14) and (3.15), we obtain 

(bll 1 - &SW(As) + (bH 1 - bi^KAs) - - ^^(A*) = . (3.16) 
Similarly, we obtain 

- biMh^As) + (&&' - ^ J Ki(A s ) - " &&')««(A*) = . (3.17) 
Therefore, from the relations (3.7), (3.8) and b\\j = b 1 ^ = 6( 7/ , we obtain 
(bi H - b[ n )a^(A s ) + (bi 11 - b{ n )a^{A s ) 



o 



xfU ;.///, 



b[ n ) \a^(A s ) + a^(As)} = 0, 



which leads to 



WW - b? 1 ) + ±b« (bi» - hi*) - 2bi\b{" - bi") 



2tt A fi 



2tt A l 



(3.18) 



= - fcf'K X (A L ) + - b?')aGr\A L ) - - &P)a£(Ai,) . (3.19) 
For the model with detMu = 0, the relation (3.19) becomes 



13a 2i? (A«) + — -ln- 



178 1 A fi 
In - — 

15 2tt A l 



127 



17 



— a3 1 (A L ) + -ai" 1 (A 



15 



6 



113 
~30 



-a 



(Aj 



(3.20) 



The right-hand side of (3.20) gives the value —97.82 for the input values oi' x (mz) 
0.01683, a L {m z ) = 0.03349 and a 3 {m z ) 

we have used the initial values at fi — mz instead of those at \x 



0.1189 [13 1, where for convenience, 

A L . The 



relation (3.20) puts a lower bound on the ratio An/ 'Al'- For a 2R (A R ) > 1, we 
obtain A R /A L > 2 x 10 135 (for A S /A R = 50 0) and A R /A L > 3 x 10 22 (for 
As/A R > 1). Such a large value of A r /Al is physically unlikely, so that the case 
A is ruled out. 

By similar discussion to the relation (3.19), it turn out that the conclusion 
that the case A is ruled out is still unchanged for the model without the condition 
detMu = and also for the minimal SUSY version of the present model. 



4 Case of SO(10) SU(2) x SU(2) x SU(4) 

Next, we investigate the case B, SO(10) L x SO(10) R -> [SU(2) x SU(2)' x SU(4)] L x 
[SU(2) x SU(2)' x SU(4)] fl . At the energy scale fj, = A s , the symmetries [SU(2)' x 
SU(4)] L x [SU(2)' x SU(4)] fl are broken into U(l) y x SU(3) C by Higgs scalars 

(1,2, 4; 1,2,4), 

(1,1, 10; 1,1,1), (4.1) 
(1,1,1; 1,1, 10), 

i n 



of [SU(2) x SU(2)' x SU(4)] L x [SU(2) x SU(2)' x SU(4)] R , where Higgs scalars 
$y, $£, and $^ generate the masses Mp, Ml and Mr, respectively. In the present 
section, we call the regions < /i < Ar, Ar < /i < As, and As < fi < Ax regions 
I, II, and III, respectively. 

The electric charge operator Q is given by Eqs. (3.2) and (3.3) in the regions 
I and II, respectively, but the relation (3.4) is replaced by 



lY = I? + ±Y L + I'* + ±Y R , (4.2) 
so that the boundary condition (3.7) is replaced by 

^\As) = ok 1 (A fl ) + l^KAs) + aft (As) + |a^(As) . (4.3) 
The boundary conditions at /i — As and /i = Ax are as follows: 

a^(As) = a^(As) + a^(A s ) , (4.4) 





= a 3 l(As) = a 4 l(A s ) , 


(4.5) 


«ri(A s ) 


= «3r( a s) = a^(A s ) , 


(4.6) 


«2l(Axl) ~- 


= ^2l( a xl) = ol^KAxl) , 


(4.7) 


®2r( A Xr) ~- 


= ®'2R\ A XR) = «4i?( A Xfi) , 


(4.8) 



where, for convenience, we have again distinguished the unification scale of SO(10) l, 
A XL , from that of SO(10) fl , A XR . 

Since b'H 1 = b'lH = b'i n ^ bH 1 = bljl = b{ n , we obtain 

a' 2 - L \A s ) - c£(As) = (&'//' - bi 1 /)^- In (4.9) 

Ztt A X l 

a' 2 - R \A s ) = a^(As) = (b'% - b^ln (4.10) 



XR 



i.e., 



a^(As) + a£(A s ) + 2{b[ n - bf 1 )^- In ^, (4.11) 



1 1 



where Ax = (AxlAxr) 1 ^ 2 ■ On the other hand, from Eqs.(4.3)-(4.6), we obtain 



« 3 _1 (A5) + 2 [^(As) + ^fl(As)] - ar^Ag) = , 



(4.12) 



so that 



« 3 '(As) + 2 [a 2 L(As) + « 2 r(A5) 



1 . A 



V (A 5 ) + 3$" - 6^ ) 7^ In ^ = 0. (4. 13) 



72 " 2 ^"A* 



Similarly, from Eq.(4.7), we obtain 



«Sl (As) " ^(As) + - 6&') ^ ^ ^ 



0, 



(4.14) 



so that, together with the equation with (L — > R) in (4.14), we obtain 



a^(A s ) - [a^(A s ) + a^(A s )} + 2(tf" - b[ H )^\n ^ = 0. (4.15) 



By eliminating A X /A R from (4.13) and (4.15), we obtain 

c 3 aJ^As) + c 2 [a^(As) + a^(A s )] - a a^(A s ) = 0, 

where 



2 -jW'-W 1 )- 



Since 



a^(A s ) 



5 i/» x 3 i.. . ,,, 1 , Aj 5,„ 1 , An 
-«r (Al) - -^(A R ) + In - + -b»- In ^, 



(4.16) 

(4.17) 
(4.18) 

(4.19) 
(4.20) 



a£(A 5 ) = a^(A L ) + 6^ In ^ + 6^ In ^, 
^(As)=« 2L 1 (A fi ) + 6^1n^, 



2tt A 



R 



(4.21) 
(4.22) 



1 o 



a^(As) = a£\A L ) + b?- In ^ + b{- In ^, (4.23) 



2tt A fi 

the relation (4.16) leads to the constraint for Ar/Al- 

= (c 2 + |d) a^(A fl ) + + 2c 2 b{ 1 - c^ 1 )^ In ^ 

+ (c 3 &3 + - ^i) ^ ln xf + c 3«3 ^l) + <%oql(k L ) - ^\A L ) 

= 19.5 a^(Afl) + 19.67 log ^ + 32.31 log ^ - 193.96 , (4.24) 

Al A r 

where we have used the values of bi given in Table 2 and the same input values 
of ai 1 (A L ), a^" 1 (A £ ,) and «g X (A L ) as those used in (3.20). For A 5 /A R = 50, the 
relation (4.24) leads to 

log ^ = 7.071 - 0.9915 aji(A fl ) , (4.25) 

so that, for a^(A K ) > 1, we obtain the constraint 

k = A R /A L < 1.20 x 10 6 . (4.26) 
Similarly, we can obtain the constraint for A x / A s : 

log ^ = 4.098 + 0.8517 aJjJ(Afl) . (4.27) 
As 

We show the values of A R , A s and Ax for the typical values of a^A^) in Table 
3. The values of A X l and A X r depend not only on the input value of a^(A fi ) but 
also on that of a^(A s ), because 

1 . A XR 



a^(A s ) = a^(As) + {b[ n - b{ n )—hi 



27T Ac 



If 1 / . 1 / . 1 / . n! /, ttt , rrr% 1 , ^-X 



2 



«m(A*) - ^l(As) + a 3 \A s ) + (b{ n - b 1 / 1 )- In 



2vr A X r 

-3.785 - 0.1964 aJ^(A fl ) + 1.405 log^- , (4.28) 



XR 
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i.e., 



log = 2.694 + 0.1398 a^(A fl ) + 0.7118 a^(A s ) , (4.29) 

A XR 

where we have used A S /A R = 50. For a^(A fi ) > 1 and a± R (As) > 1, the relation 
(4.29) gives the constraint 

Axl/Axr > 1-26 x 10 7 . (4.30) 

The relation (4.29) concludes that a model with Axl = Axr is ruled out. Values 
of A X r and A X l for typical values of a^(A^) and a~[ R (A s ) are also listed in Table 
3. 

Considering the present results |14j of the experimental searches for the 
right-handed weak bosons, we take k = A R /A L > 10, so that we conclude that 
the allowed regions of k, the intermediate energy scale As and the unification scale 
Ax = (A XL A XL ) l l 2 are 

k= 1.3 x 10 1 - 1.2 x 10 6 , 

A 5 = (6.0 x 10 4 - 5.5 x 10 9 ) GeV , (4.31) 

A x = (9.8 x 10 13 - 4.9 x 10 14 ) GeV , 

corresponding to the values a^(A fi ) = 6—1. Behaviors of the gauge coupling 
constants for a typical case are illustrated in Fig. 1. 

5 Conclusions 

In conclusion, in order to examine the idea that the extra fermions F R and Fl 
in the universal seesaw mass matrix model, together with the conventional three 
family quarks and leptons fi and f R , are assigned to (fi + F^) ~ (16, 1) and 
(fa + Fl) ~ (1, 16) of SO(10) L x SO(10) fl , we have investigated the evolution of 
the gauge coupling constants and intermediate mass scales. The case A, SO(10) L x 
SO(10) fl — > [SU(5) x U(1)']l x [SU(5) x \J(1)'] R , is ruled out because the results are 
inconsistent with the observed values of the gauge coupling constants at /J, — mz- 
The case B, SO(10) L x SO(10) B -> [SU(2) x SU(2)' x SU(4)] L x [SU(2) x SU(2)' x 
SU(4)]ij, is allowed for the intermediate energy scale A R ~ (10 1 — 10 6 ) GeV if we 
accept a model with Axl 7^ A XR , where Axl and A XR are the unification scales 
of SO(10) L and 80(10)^, respectively. We have obtained the allowed regions k ~ 
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K^-IO 6 , A s ~ (6xl0 4 -6xl0 9 ) GeV, and A x = (A X l^xr) 1/2 = (5 x 10 14 - 10 14 ) 
GeV correspondingly to a2 R (An) ~ 6 — 1. 

In the case B, since Ml ~ Mr ~ M x ~ Mp (F ^ N), the case gives effective 
neutrino mass matrix M(j/l) — —nT'LM£ 1 mJ / , so that the conventional neutrino 
masses rn(z/j) are of the order of m{ej)J ' k. However, for the condition a^(A^) > 
1, which is a condition that the model is perturbative, the value of k has been 
constrained by (4.26), i.e., k < 1.20 x 10 6 . This suggests that m{y T ) ~ m(r)//t > 
10 3 eV. Such a large value of m(v T ) is unlikely. Therefore, the straightforward 
application of the case B to the neutrino mass generation scenario is ruled out. 

However, the numerical results in Sec. 4 should not taken rigidly, because 
the calculation was done at one-loop. Moreover, the results are dependent on the 
input value Ar/ 'As- The value Ar/ As = 0.02 have been quoted from Ref. 
where the value was determined from the observed value of m c /m t on the basis 
of a specific model for ttil, ttlr and Mp. Exactly speaking, the value 0.02 means 
yLVLljRVR/ysVs = 0.02, where y's and f's are the Yukawa coupling constants and 
vacuum expectation values, respectively. Because of the numerical uncertainty of 
i/L, Vr, and ys, the numerical results may be changed by one or two order. The 
case B cannot still be ruled out. 

In the present paper, the cases for SUSY version of the model have not been 
investigated systematically, because many versions for the energy scale of the SUSY 
partners of the super heavy fermions F can be considered. Nevertheless, the case 
A can easily be ruled out by simple consideration. On the other hand, for the case 
B, it is a future task whether the SUSY version is allowed or not. 

When we take the numerical result of the constraint (4.26), we can consider 
a minimum modification of the case B. In the case B, the Dirac mass matrix Mp, is 
generated by the Higgs scalar ~ (1, 2, 4; 1, 2, 4) of [SU(2) x SU(2)' x SU(4)] L x 
[SU(2) x SU(2)' x SU(4)]a, while the Majorana mass matrices Mi and Mr are 
generated by the Higgs scalars $l ~ (1, 1, 10; 1, 1, 1) and &r ~ (1, 1, 1; 1, 1, 10), 
respectively. We assume that the symmetries SU(4) L and SU(4) fl are broken into 
[SU(3) x U(1)] L and [SU(3) x U(1)] R at y, = A NL = 0(M L ) and y = A NR = 
0(M R ), respectively, and the energy scales A NL and A NR sufficiently larger than 
As = 0(Mp>), at which all the fermions F (not /) have Dirac masses Mp and the 
symmetries SU(3) L xSU(3) R and U(l) L xU(l) R are broken into SU(3) L+/? = SU(3) C 
and \](1) L+R = U(l) y , respectively. Then, the neutrino mass generation scenario 
is changed from the scenario (b) to the scenario (a). Although the expression of 
M v is still given by M u ~ — m^Mj rrt^, the suppression factor for neutrino masses 
is changed from 1/k to (1/ k)(A s / A NL ). By taking A5/A7VL ~ 10~ 3 , we can obtain 
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reasonable values of the neutrino masses for the case a^(A^) ~ 1. Of course, in 
the modified version with A X l A NL ^> A s , the unification scales of A X l and 
Axr are changed by an order of one or two. However, Ar and As are insensitive 
to the present modification. 

In the present paper, we have not discussed the evolution of the Yukawa cou- 
pling constants. The phenomenological success in Ref. |I| has been obtained by tak- 
ing b e = 0,b u = -1/3 and b d = -e iPd (f3 d = 18°), where M F = m A/diag(l, 1, 1 + 
36/) in the basis on which Mp is diagonal. The shapes (not the magnitudes) of 
Me = moA e diag(l, 1, 1) and Mu = moA u diag(l, 1, 0) are almost invariant under 
the evolution, while the shape of ~ m Arfdiag(l, 1, —2) is not invariant. The 
following problems remain as our future tasks: (i) what value of bd is favorable at 
the unification scale fi = Ax', (ii) whether we can still assert X u ~ A^ or not; (iii) 
whether the mass matrix itlr can still be approximately diagonal on the basis on 
which is diagonal; and so on. The numerical results in Ref.[l[] will be somewhat 
cahnged under the present SO(10) L x SO(10) R model. 

In any case, for the universal seesaw mass matrix model based on the SO(10) L x 
SOflO)^ unification, if we consider the symmetry breaking SO(10) L x SO(10) fi — > 
[SU(2) x SU(2)' x SU(4)] L x [SU(2) x SU(2)' x SU(4)] fl and we accept the case 
Axl 7^ where A X l and A R are the unification scales of SO(10) L and SO(10) R , 
respectively, we can find a solution of the intermediate energy scales and As 
for the unified description of the quark and lepton mass matrices, where only the 
top quark mass mt is given by mt ~ Al in contrast with n g < (g/t) and the 
neutrino masses m(i>i) are reasonably suppressed compared with the charged lep- 
ton masses m(ej). The model is worth while being taken seriously as a promising 
unified model of the quarks and leptons. 
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Table 1: Quantum numbers of the fermions / and F and Higgs scalars 4>l, 4>r 
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Table 2: Coefficients in the evolution equations of gauge coupling constants. The 
cases A and B are cases with the symmetry breaking patterns SO(10) — > SU(5) x 
U(l) and SO(10) -> SU(2) x SU(2) x SU(4), which are discussed in Sees. 3 and 4, 
respectively. 
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| ^ = -53/10 
| b 1 " = -53/10 


| 6' 2L = -13/6 
| b' 2R = -13/6 



Table 3: Intermediate mass scales A R and A s versus a 2R ~(A R ) in the case of 
SO(10) -> SU(2) x SU(2) x SU(4). Input values A R /A S = 0.02 and A L = m z = 91.2 
GeV are used. The upper and lower rows of A X r and A X l correspond to the values 
for a^ R (A s ) = 1 and a^ R (A s ) = 2, respectively. 
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Figure 1: Behaviors of (dotted line) in A L < /i < A^, af^/i) (dotted 

line) in A R < /i < A s , a^KlA (solid line) in A L < fi < A X l, o^mif 1 ) (solid line) 
in A R < /i < Axr, «3 x (/i) (dashed line) in A^ < /i < As, o^l 1 ^) (dotted line) in 
As < /i < Axl, and Ql I 2r{i j ) (dotted line) in As < /i < Ax_r, aj^/i) (dotted chain 
line) in As < /i < Axl, and a^ R (fi) (dotted chain line) in As < /i < Ax_r, where 
A L = 91.2 GeV, A R = 1.10 x 10 8 GeV, A s = 5.48 x 10 9 GeV, A XR = 1.39 x 10 11 
GeV and A XL = 1.71 x 10 18 GeV. The values c/f : (A L ) = 59.42, ^(Al) = 29.86, 
« 3 _1 (A i ) = 8.410, a^(A^) = 1 and a^ R (A s ) = 1 are used as the input values. 
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